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ABSTRACT 

We discuss the deflection of light and Shapiro delay under the influence of gravity as described by Schwarzschild 
metric. We obtain an exact expression based on the coordinate velocity, as first set forth by Einstein, and present 
a discussion on the effect of velocity anisotropy. We conclude that the anisotropy in the coordinate velocity, 
as the velocity apparent to a distant observer, gives rise to a third order error in the deflection angle, so that 
the practical astronomical observations from gravitational lensing data remain inconclusive on the anisotropy. 
However, measurement of Shapiro delay provides a fairly convenient way to determine whether the spacetime is 
optically anisotropic for a distant observer or not. We calculate the Shapiro delay for a round trip path between 
Earth and Venus and observe excellent agreement to two experimentally reported values measured during a time 
span of six months in 1967, without any need to extra fitting parameters. This is while the expected delay 
obtained from an isotropic light velocity as described by Einstein's model suffers from much larger errors under 
similar conditions. This article illustrates the usefulness of the equivalent medium theory in understanding of 
general theory of relativity. 
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1. INTRODUCTION 

It is widely known that Schwarzschild metric is the only exact solution of the Einstein field equations under 
spherical symmetry, which is given by the metricEB 

cldr 2 = c 2 (l - ^)dt 2 - - r 2 (d0 2 + sin 2 <j>d9 2 ), (1) 

r 

where r is the proper time, Co is the speed of light in flat space, (r, 6, cj>) constitute the spherical polar coordinates, 
and r s is the Schwarzschild radius of the star. Initially in 1911, Einstein mentioned that the light velocity arising 
from such a spherically symmetric gravity should be modified aJ3 



c(r) = co [1 + *(r)] , (2) 

in which $(r) = — ^ is the gravitational potential. He furthermore took notice of this radial dependence of light 
velocity to propose the possibility of deflection of light while passing across a massive object according to the 
angle 

(3) 

with R being the closest distance of approach of the light ray with respect to the center of the massive object. 
But later in 1916, Einstein noticed a missing correction factor of 2 in ^ as 

c(r) = c [1 + 2$(r)] , (4) 
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resulting in 



Famously, this result was confirmed in the 1919 experiment of Sir Arthur S. EddingtonP 

Perhaps, Eddington was the first to point out explicitly that "For the original coordinates of ([I]) the velocity 
of light is not the same for the radial and transverse directions" P To resolve this difficulty, he proposed the 
change of coordinates using the transformation 



r=(l + ^) 2 p, (6) 
which takes advantage of the non-physical radial coordinate p, and subsequently results in the isotropic metric 



c o^ 2 = - (1 + W + p 2 # 2 + p 2 sin 2 0^ 2 ). (7) 

It is readily possible to calculate the resulting coordinate velocity of light in this frame to be 

1 _ ti 

c(r) = c (i + | ? )3 ^c [1 + 2$(r)] , (8) 

where the approximation r ~ p holds in the limit of p >> r s , according to ([6]). The light velocity Q as measured 
from a distant observer is apparently isotropic, while a local observer regardless of the choice of coordinates would 
always see a constant and direction- independent speed of cq , which is of course also independent of the choice 
of coordinates. This also follows the principle of equivalency, which states that any accelerated observer would 
always see a flat space to a first-order of accuracy. 

As a matter of fact, it is easy to observe that there is a maximum anisotropy in the coordinate speed of light 
as suggested by ([!]) equal to |$(i?)|co, corresponding to the difference between radial and transverse speeds, the 
exact relation of which is found through the equivalent medium theory to bd33 



c(r,k) = c a/1 + 2$(r) a/I + 2$(r) cos 2 #(r,k). (9) 

Here, ^(i^k) = arccos(f ■ k) is the angle made by the unit radial f = — and unit propagation vectors k = ^, and 
hence is anisotropic as seen by a distant observer. Interestingly, in the limit of large radius r, ([9]) simplifies to 



c(r,k) « c {l + [l + cos 2 (r,k)] $(r)}, (10) 

which establishes the fact that Einstein's correction coefficient of 2 indeed varies between 1 and 2 depending 
on the angle of propagation. Evidently, this value remains close to 2 as long as the light ray is at a position 
r away from its closest distance of approach R; it only for the region sufficiently close to the closest distance 
of approach, appreciably decreases to achieve a minimum of 1, exactly at r = R. The above result |10j), which 
is obtained directly through the equivalent medium theorypH has been also deduced by Eckstein^ through 
geometric considerations. It is furthermore noticed that the corresponding anisotropy in the light velocity due 
to the gravitational field of sun and the black hole at the center of milky way galaxy could respectively reach 
9 x lO- 9 P0and 5 x 1CT 7 according to MizushimaP 

In 1959, Yilmaz^l proposed an experiment which put the equivalence principle of the general relativity to 
the test. He suggested to set up an optical interferometer, looking for any such possible local anisotropy arising 
from the gravity of Sun. Ten years later, Shamir and Fox published a papei^ reporting null results within an 
accuracy of 3 x 10~ n . 

The purpose of this manuscript is two fold. We examine the deflection of light as well as Shapiro delay as 
plausible methods to detect the possible anisotropy in the light speed: 



• First, we show that the exact deflection angle arising from an anisotropic light velocity (|9| could mean- 
ingfully differ from the value as predicted from the isotropic velocity Q. For this purpose, we first derive 
an exact expression for the deflection angle as suggested from Q. This expression agrees with the more 
accurate expression reported recently by VirbhadrElHJli^l to the second-order. Then we numerically evaluate 
the exact angle of deflection based on ^ to notice an appreciable difference when R < 5r s . We discuss 
that for all practical reasons the deflection of light is not a good means to decide on the anisotropy of 
equivalent medium space, since the expected difference is vanishingly small when R » r s . 

• Secondly, we consider the Shapiro delay. We calculate this delay based on anisotropic description of space 
and observe an excellent agreement to the experimental data, without any need to an external parameter. 
This is while either of the 1911 or 1916 Einstein's isotropic expressions suffer from much larger errors. 



2. DEFLECTION OF LIGHT 
2.1 Isotropic Coordinate Velocity 

The significant difference between the concepts of coordinate velocity and local velocity of light has been a matter 
of debate in numerous works, which is summarized adequately by LonghiJ-^ More recently, it has been stated by 
JenserP^ that the Schwarzschild metric could violate the weak principle of equivalence, which is here understood 
to be a misinterpretation of the concept of coordinate velocity. Similarly, it has been arg ued by the authoJZEU 
that whether local interferometry could ultimately detect gravitational waves or not; this is here left as the 
subject of a separate study. 

In order to calculate the angle of deflection, we follow the integral method originally used by Einstein.E For a 
light ray moving on the plane (xi,X3) subject to the gravitational field of a massive object located at the origin, 
the deflection angle a can be found as^ 

OL= \ —. r — dx 3 , 11 

J-oo c(r) dxi 

with r = (xi, X2, X3) and X2 = 0. Upon plugging ^ in (11) we get after significant algebra 

7r — ArctanV a 2 — 1 tt 

a{a) = a -, (12) 

ya — 1 * 

in which the definition a = is adopted. In the limit of R — > 2r s the angle a approaches infinity, which can 
be interpreted as infinite rotation at R = 2r s . Since the photon sphere is known to be located at R = |r s , the 
expression ( 12 ) might be regarded as invalid for a close to unity. This behavior is illustrated in Fig. [lj 



When a > 1, this expression may be expanded in the powers of - to obtain 



1 tt 1 2 1 3tt 1 8 1 
a(a) « - + -— + -— + — — + — — . (13 
a A a 2 3 a 3 16 a 4 15 a 5 



This agrees to the Einstein's result |5| to the first order. In contrast, Virbhadra's expression read^MH 

1 1 , 157T ,1 . 

a 4 16 a z 

which has been obtained from direct integration of an extended expression by Weinberg^ to the case of the most 
general static spherically symmetric asymptotically flat universe. 



Angle of Deflection (in unit of 2n) 




2.2 Anisotropic Coordinate Velocity 

Now, we proceed to calculate the deflection angle as caused by using the exact anisotropic coordinate velocity 
|9). We may use 

+ °° 1 (15) 



c(r,k) dx\ 

The difficulty in evaluating the integral is that the propagation vector k changes its direction along the light 
trajectory. Hence, we first rewrite (15 1 as 

°«* 2 f^ a tr w ,16) 

Upon differentiating with respect to X3 from both sides we arrive at the first order nonlinear differential equation 

^ ( ) ~ 2 ^ ^ C i r i^ L ) 

dx3 c(r,k) Qxi 

in which we have 



fc = cosa(r)ii + sina(r)i3. (18) 

Now, using the initial condition a(R, 0, 0) = we can obtain the deflection angle as a — lim z _ i . 00 a{R, 0, z), 
which can be easily done by means of accurate numerical integration methods. 

2.3 Comparison 

The variation of deflection angle a versus the closest distance of approach R is here plotted as shown in Fig. J2J 
based on the four expressions. Dotted, dashed, and dot-dashed curves respectively represent the Einstein (15]), 



Virbhadra (14 1, and our expression (12). In contrast, the solid heavy curve represents the angle of deflection 
( 15 ) based on anisotropic coordinate velocity. 

In general, the difference between the deflection angles from anisotropic and isotropic coordinate velocities 
is less than 1% when a > 5. So, observation of angles of deflection may not be an efficient way to determine 
whether the space behaves optically anisotropic or not. 
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Figure 2. Angle of deflection from various expressions versus normalized parameter a. 



3. SHAPIRO DELAY 

Calculation of Shapiro delajU^MD j g a sur p r i s i n gly convenient method, as we here show, to determine whether the 
space looks optically isotropic for a distant observer or not. In 1964, Irwin I. Shapiro proposed 18 an interesting 
measurement to validate the general theory of relativity, through observation of round trip delay of a radio signal 
traveling between Earth and Venus, while the radio beam grazes the solar limb. Since all the existing relations 
for light velocity are independent of frequency, the same equations can be applied to radio waves. Under the 
influence of Sun's gravity, the speed of the electromagnetic wave would decrease. As a result of this slowness, a 
time-difference can be detected by subtracting the round trip delay of a radio beam propagating at the constant 
speed of Co, from the actual measured value, which can be attributed to the influence of the gravitational field. 
Based on an isotropic light velocity as suggested by Einstein in 1916 Q, Shapiro's predicted value turns out to 
be about 200/xsecJ^ while using the Einstein's 1911 result ^ one may expect half as much, or roughly 100/zsec. 
Soon after, Ross and SchifP- noticed the dependence of calculations on the particular choice of coordinates of 
either ([T]) or ([9]), and suggested the correct transformation to eliminate the apparent coordinate dependence. 
More recently, Shapiro delay has been used to accurately measure the mass of a Neutron starP2l 

Two measurements actually took place in May and September 1967,^ respectively giving rise to the values of 
160/isec and 150/isec, as illustrated in Fig. [3j Shapiro's observations^! matched the theory within the relatively 
large error of ±20%, while using a generalized metric with |(1 + 7) = 0.9 ± 0.2, instead of the Schwarzschild's 
with 5(1 + 7) = 1- It has been more recently shown through later astronomical observations of gravitational 
lensing^SJ by Cassini spacecraft, that this expression ^(1 +7) should be very close to unity, indeed within an 
error margin of only ±10 -5 . 

It is fairly easy to make an estimate of the Shapiro delay based on the anisotropic coordinate velocity 
We may here write down 



At: 





' 1 1 


1 —ClT 


_c c(r,k)_ 



dx 3 , 



(19) 



in which ot and an are respectively the orbiting radii of Earth and Venus around the Sun, being equal to 
1.496 x 10 n m and 1.0829 x 10 n m. The closest distance of approach R is taken to be the solar radius 6.96 x 10 8 to 
having a Schwarzschild radius of approximately 3 x 10 3 m. Numerical evaluation of the integral ( 19 ) gives rise 



to the value of 148.342/isec, in surprisingly good agreement with the observed values of 150 — 160/isec. This is 



200 




while using the 1911 ^ and 1916 Q expressions give respectively the values of 118.113/isec and 236.226^sec. It 
should be added here that in this calculation, we have assumed circular orbits for Venus and Earth, while their 
actual orbits have eccentricities of respectively 0.68% and 1.7%. This agreement also narrows the experimentally 
determined value of 7 = 1 through radio Shapiro delay measurement to less than 1%, thus reaching to the best of 
similar observations so farPil It is noticeable that EcksteirPhas also used (10) to obtain a confirming expression 
to that of the generally accepted formula inPS 

This may provide a very good basis for the understanding of how coordinate and local velocities of light 
differ in practice. Actually, the coordinate velocity becomes meaningful in the framework of equivalent medium 
theory, and we have shown the accuracy of this approach as an alternative calculation method. The calculations 
take place on a standard Minkowskian or Cartesian frame of coordinates, instead of the curved space. Moreover, 
this reasonable agreement between our calculation and observed data reveals the good accuracy of the general 
theory of relativity, without any need for any externally inserted fitting parameters. However, this observation 
relies on the fact that the velocity of light as appears to a distant observer is actually anisotropic under the 
influence of gravitational field. In short, we may deduce that the coordinate velocity, or the apparent velocity to 
a distant observer who is outside a gravitational field, actually corresponds to the same quantity obtained from 
the standard equivalent medium theory. 



4. CONCLUSIONS 

In this paper, we discussed the possible influence of gravity-induced anisotropy on the propagation of optical 
waves within the framework of equivalent medium theory. We discussed two plausible approaches to examine the 
physical influence of this anisotropy through measurement of deflection of light and also the so-called Shapiro 
delay. We concluded that for practical reasons, the deflection of light is not sufficiently accurate to distinguish 
between isotropic and anisotropic coordinate velocities of light in the equivalent medium framework. On the other 
hand, measurement of the Shapiro delay turns out to be an ideal method to test the anisotropy of coordinate 
light velocity, as it appears to a distant observer. We found a fairly good agreement between the calculated value 
based on an anisotropic light velocity and the experimentally observed value, being within 1% and 7% of the 
observed values, belonging to two separate measurements. 
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